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Abstract 

We propose a new type of reduction for integrable systems of coupled matrix 
PDEs; this reduction equates one matrix variable with the transposition of another 
multiplied by an antisymmetric constant matrix. Via this reduction, we obtain a 
new integrable system of coupled derivative mKdV equations and a new integrable 
variant of the massive Thirring model, in addition to the already known systems. We 
also discuss integrable semi-discretizations of the obtained systems and present new 
soliton solutions to both continuous and semi-discrete systems. As a by-product, a 
new integrable semi-discretization of the Manakov model (self-focusing vector NLS 
equation) is obtained. 
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1 Introduction 



Since the seminal work of Manakov pQ in the early 70s, integrable systems of coupled 
partial differential equations (PDEs) that are associated with higher than second-order 
matrix spectral problems (Lax pairs) have been the focus of intensive research. In par- 
ticular, of prime importance among such systems are the vector PDEs that are invariant 
under the action of a classical matrix group on the vector dependent variables; this invari- 
ance also represents a "symmetry" or gauge invariance of the spectral problem. Because 
of this large "symmetry", the vector PDEs usually allow the existence of solitons with 
internal degrees of freedom that exhibit highly nontrivial and interesting behaviors in 
soliton interactions. Moreover, their simple and symmetric form of equations very often 
leads to their potential or practical applicability in various branches of physics as well as 
applied mathematics. Typical examples of such vector PDEs include the U(m) -invariant 
Manakov model [HE], also referred to as the self-focusing vector nonlinear Schrodinger 
(NLS) equation, 

m 

i^t + <3xx + 2||g|| 2 g = 0, q = (qi,q 2 ,...,q m ), \\q\\ 2 := q ■ q ] = ^ \qj\ 2 ; (1.1) 

i=i 

two distinct versions of the vector mKdV equation [3H5] with 0(m)-invariance, 

m 

qt + Qxxx + 6(q,q)qx = Q, <?= 0i,?2,---,gw), (q,q) ■= q- q T = ^2q], (1.2) 

3=1 

q t + q xxx + 3(9, q x )q + 3(q, q)q x = 0; 
the vector third-order Heisenberg ferromagnet model with 0(m)-invariance [6j[7], 

S t + S xxx + ^((S x , S X )S) X = 0, (S, S) = 1; (1.3) 

and the 0(m)-invariant vector extension [7] of the third-order Wadati-Konno-Ichikawa 
equation [8], 

r q x 



q t + 



o. 



The recent developments in computer algebra packages and improvements in CPU per- 
formances have further increased the number of m-component integrable systems with 
U{m)- or 0(m)-invariance as well as (m + l)-component systems with 0(m)-invariance 
to a considerable extent [9HTT]. However, in contrast to the U(m) -invariant and 0(m)- 
invariant systems, little research has been conducted on integrable systems having invari- 
ance with respect to the symplectic group Sp(m); this refers to the group of 2m x 2m 
real/complex symplectic matrices, Sp(m, M) or Sp(m, C), in accordance with the attribute 
of the dependent variables. To the best of the author's knowledge, the only example of 
a 2m-component vector nonlinear PDE with Sp(m) invariance is the system of coupled 
derivative mKdV equations studied using the bilinear method by Iwao and Hirota [12J, 



du i d^Ui 
dt dx 3 



III r\ n 

E/0U 2 j-l OU 2 j 

\ dx 2j 2j ~ dx 

.3=1 



^=0, z = l,2,. ..,2m, (1.4) 
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as well as its higher symmetries. It should be noted that system (jl.4p is a natural multi- 
component generalization of the two-component system [ (jl.4j) with m = 1] derived within 
the framework of the Sato theory by Loris and Willox [121 E]. Here, the Sp(m) in- 
variance refers to the fact that system (ll.4p is form-invariant under the following linear 
transformation: (ui, u 2 , . . . , u 2m ) l— {ui, u 2 , ■ ■ ■ , u 2m ) S T , where S is an (x, /^-independent 
element of the symplectic group Sp(m) with a proper ordering of the base vectors, and 
the superscript T denotes the transposition. 

The main objective of this paper is to expand the class of 5p(m)-invariant integrable 
systems and to characterize some of their interesting properties within the framework 
of the inverse scattering method. To achieve this goal, it must be first noted that sys- 
tem (II. 4p is homogeneous with respect to the following weighting scheme: w(d x ) = 1, 
w(d t ) = 3, and w(ui) = 1/2; this weighting scheme is the same as that for the third-order 
symmetries of vector derivative NLS (DNLS)-type systems. Motivated by this observa- 
tion, we start with more general systems, that is, the third-order symmetries of integrable 
matrix generalizations of the DNLS-type equations. Then, we propose a new type of re- 
duction for these third-order integrable matrix PDEs, at least in its explicit form, wherein 
one matrix variable is related to the transposition of the other multiplied by an antisym- 
metric constant matrix [T5]. Considering the special case wherein the matrix variables 
are restricted to the form of column/row vectors, we obtain two Sp(m)-invariant systems; 
one coincides with system ( II. 4p . while the other one, 



dui d 3 Ui d 
dt dx 3 dx 



E 



du 



2j-l 



du 



dx 



u 2j 



u 2j - r 



2/ 



dx 



U; 



0. 



1,2, 



2m, 



(1.5) 



appears to be a new integrable system [15] . Note that the location of d/dx in (I1.5P is 
different from that in (jl.4p . 

Once we have identified systems (II .4p and ( \1.5h as the reductions of the matrix DNLS- 
type systems, it is not difficult to further extend the class of Sp(m)-invariant integrable 
systems. First, we consider the integrable matrix generalizations [TBI ITT] of massive 
Thirring-type models [181422] . which are the first negative flows of the matrix DNLS hi- 
erarchies. Then, via the same type of reduction, we can directly obtain a new integrable 
variant of the massive Thirring model that is a hyperbolic system with 5p(m)-invariance, 



d 2 Vj 
drdx 



+ Vi 



E 



dv' 



2j-l 



dx 



v 2j 



v 2j -r 



dv' 



2.7 



dx 



0. 



1,2, 



, 2m, 



1.6) 



as well as an equivalent system up to the interchange of x and r. Second, we consider 
the integrable space discretizations (semi-discretizations, for short) of the matrix DNLS 
hierarchies, including the massive Thirring-type models, proposed in ref. [23] Though not 
all of the semi-discrete matrix DNLS-type systems in ref. [23] are useful for our purpose, 
we find that the semi-discrete Kaup-Newell hierarchy [cf. (3.1) and (5.2) in ref. [23] allows 
proper reductions to yield integrable semi-discretizations of the above 5*p(m)-invariant 
systems. Note that the discrete analogue of the reduction used to obtain a continu- 
ous 5*p(m)-invariant system may not be uniquely determined. In fact, in addition to 
a single integrable semi- discretization of system (II. 4p . we obtain two integrable semi- 
discretizations for each of the Sp(m)-invariant systems (II. 5p and (jl.6p . This result partly 
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illustrates the wide applicability of the semi-discrete Kaup-Newell hierarchy provided in 
ref. |23]to the theory of integrable discretizations. 

Hereafter, we will not restrict ourselves to the canonical representation of the Sp(m)- 
invariant systems, but will present them in a slightly generalized (but still being inte- 
grable) form. Specifically, we replace terms such as Y^=i{ ~^faT~ u 2j ~ u 2j-i^r) by those 

such as Ek^xkm^)'^^^ — u olhf)i wnere the integer M may or may not be even, 
and Cjk (j < k) are arbitrary coupling constants (cf. ref. H2]). It can be recalled that 
according to the classical theory of matrices (see, e.g., refs. I2HI25]), any antisymmetric 
matrix C := (Cjk)j,k=i,...,M can be transformed to the block diagonal form 




J 



1 
-1 



with an invertible matrix P. Here, 2m (< M) is equal to the rank of C , where m denotes 
the number of copies of J. This guarantees that through an invertible linear change of the 
dependent variables, the generalized system involving Cjk can always be converted into 
an Sp(m)-invariant (sub)system in the canonical form with linear equation(s) coupled to 
it, if any (the case of 2m < M). Therefore, in this paper, we also use the term "symplectic 
invariance" for such generalized systems. 

Using an approach based on the inverse scattering method, we can construct multi- 
soliton solutions of the matrix DNLS hierarchies under appropriate boundary conditions 
in both the continuous and discrete cases [26] . Thus, the soliton solutions for the Sp(m)- 
invariant systems contained in these hierarchies can be obtained by properly specializing 
the soliton parameters involved in the solutions for the latter. Moreover, as one of the 
main advantages of the inverse scattering method, the most general soliton solutions are 
obtained under the specified boundary conditions. However, the computations and discus- 
sions required to arrive at simple explicit formulas for the multi-soliton solutions are rather 
extensive and involved; hence, we do not present them here. We will defer the detailed 
derivation and investigation of the multi-soliton solutions to a subsequent publication [26] . 
In this paper, we assume decaying boundary conditions at spatial infinity and start with 
the linear integral/summation equations associated with the continuous/discrete matrix 
DNLS hierarchies, omitting their derivation via the inverse scattering method. These lin- 
ear integral/summation equations, referred to as the Gel'fand-Levitan-Marchenko type, 
provide an exact linearization [27\ of each nonlinear system under study in that they 
provide a relation between the solutions of the nonlinear system and those of the corre- 
sponding linear system. Considering a special case with a proper reduction, we solve the 
Gel'fand-Levitan-Marchenko equations to obtain the one-soliton solutions of the Sp(m)- 
invariant systems under the decaying boundary conditions. It should be emphasized that 
the soliton solutions of (11. 4p obtained in this manner are indeed more general than the 
previously known solutions [12],[13] . Although the accuracy of these one-soliton solutions 
can easily be verified by direct substitutions, unlike the case of the NLS equation, it 
is not easy to obtain such solutions directly without resorting to the inverse scattering 
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method or other sophisticated methods in soliton theory. Indeed, any naive ansatz for 
the travelling wave solutions, e.g., a complex plane wave modulated by a real envelope 
moving with a constant velocity, is most likely to fall into a trivial subclass of the most 
general one-soliton solutions, such as "the one-soliton solution" proposed by Loris and 



This paper is organized as follows. In section 2, we demonstrate that the continuous 
systems ( 11.41) . (11.51) . and ( 11.61) in a slightly generalized form, as mentioned above, are 
obtained through a new type of reduction of the matrix DNLS hierarchies. We also 
present their bright one-soliton solutions. In section 3, we propose the integrable semi- 
discretizations of these continuous systems and present the one-soliton solutions for the 
most interesting semi-discrete systems. We also derive a new integrable semi-discretization 
of the Manakov model (II. ip from one of these semi-discrete systems. The last section, 
section 4, is devoted to concluding remarks, wherein we state that the reduction considered 
in this paper is not restricted to the DNLS-type systems, but is also applicable to a matrix 
generalization of the Yajima-Oikawa hierarchy. 

2 Reductions of continuous matrix DNLS hierarchies 

In this section, we propose the Sp(m)-invariant integrable systems via reductions of the 
continuous matrix derivative NLS hierarchies. We also present the Lax pairs, associated 
linear integral equations, and one-soliton solutions for the obtained systems. 

2.1 Coupled derivative mKdV equations of Chen— Lee— Liu type 

2.1.1 Lax pair for the second flow of the matrix Chen— Lee— Liu hierarchy 

A derivative nonlinear Schrodinger (DNLS) equation iq t2 + q xx ± ilgpq'z = 0, which is often 
referred to as the Chen-Lee-Liu equation [28], permits an integrable matrix generaliza- 



Here, q and r are l\ x l 2 and l 2 x Z x matrices, respectively. Note that O on the right-hand 
side of the equations implies that the dependent variables can take their values in matrices. 
The lower indices of t are used to distinguish between different (and commutative) time 
evolutions; however, in the following text, we often omit these indices for brevity. In 
this paper, we are more concerned with the next higher flow in the matrix Chen-Lee-Liu 
hierarchy that commutes with the first nontrivial flow (12.11) . It is written as (up to a 
scaling of the time variable) [TBII29, 30j 



Willox [13J. 



tion HSKTHESIEDI 




(2.1) 




(2.2) 
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The Lax pair for (I2.2p is given by 



u = K 2 
v = K" 

+ iC 

+ c 

+ i 



■h 



-AI X 



+ c 


Q 




" 


+ i 




r 




\ r( l . 



4Jo 



c 5 



iq 



4r 



+ K 4 



-2qr 



2rq 



+ c 3 



(2.3a) 
2iq x + qrq 



-2ir x + rqr 



-i(q x r - qr x ) - \{qr) 2 



\(rq x - r x q) + \{rq) 2 



-Qxx + HVxrq - qr x q + 2qrq x ) + j(qr) 2 q 



-r xx - \Qr x qr - rq x r + rqr x ) + j(rq) 2 r 







-k(r xx q - r x q x + rq xx ) - -{2r x qrq - rq x rq + rqr x q - 2rqrq x ) + \irqf 



(2.3b) 



Here, ( is the spectral parameter independent of x and t; I\ and 1% are the l\ x l\ 
and I2 x I2 unit matrices, respectively. Let us substitute the Lax pair (12. 3 p in the zero- 
curvature condition 



U t - V x + UV - VU = O, (2.4) 

which is the compatibility condition for the overdetermined system of linear PDEs, 

# x = UV, % = VV. (2.5) 

Subsequently, equating the terms with the same powers of ( to zero, we obtain the third- 
order matrix Chen-Lee-Liu system (12. 2p without any contradiction or additional con- 
straint. 



2.1.2 Reduction 

Both the first flow (12.11) and the second flow (12.21) of the matrix Chen-Lee-Liu hierarchy 
permit the reduction of the Hermitian conjugation [33] r = Aiq^A2, where A\ and A2 
are constant Hermitian matrices: A] = A and At + = A; r = O. Moreover, unlike the first 



flow (12. ip . the second flow (12.21) allows an interesting reduction such that r is identically 
equal to Cq T , where C is an antisymmetric constant matrix: C T = —C, Ct = C x = O. 
In fact, the reduction r = q T C (or, more generally, r = Bq T C, B T = B, C T = —C) can 
also be considered for (12.21) . but we do not exploit this reduction in order to maintain a 
natural and easy-to-read flow of the paper. For the reduced system to assume a concise 
form without the imaginary unit and fractions, we consider the following vector reduction 
(h = l,l 2 = M): 



q = (ui, . . . , u M ), r = 2\Cq T = 2i 



Sfe=l CMkUk J 



(2.6) 
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which also implies the simple relation qr = 0. System (I2.2p is then reduced to a system 
of coupled derivative mKdV equations [T2] : 



du,; d 3 



U: 



dt ^ dx 3 



+ 3 



J2 a fdUj 



dx 



U k - Uj 



duk 
dx 



duj 
dx 



0. 



z = l,2, 



M. 



(2.7) 



_l<j<k<M 

Using the vector notation, (12.71) can also be written as 

u t + u xxx + 3(u x C, u)u x = 0, C T 



-C. 



Here, u = (ui, U2, ■ ■ ■ ,um) and {u x C, u) = u x Cu T . The Lax pair for this reduced system 
is obtained by substituting (12.61) into q and r in (12. 3p ; using a simple gauge transformation, 
the Lax pair can be rewritten in the form 



u = 


—A u 
\Cu T -Cu T u \ ' 




(2.8a) 




\ 3 + 2\(u x C,u) 


—\ 2 u + \u x 


- u xx - 2(u x C, u)u 


V = 


-X 3 Cu T - \ 2 Cu T x 

- \Cu T xx - 2X(u x C, u)Cu T 


X 2 Cu T u + \C(ulu 
+ C(ul x u - u T x u x ^ 


- u T u x ) 

- u T u xx ) + 2(u x C, u)Cu T u 



(2.8b) 

where A := 2i( 2 is the "new" spectral parameter. In the simplest nontrivial case of M = 2, 
setting 

C12 = i, ui = ip, u 2 = ip*, 
we obtain the following single equation [T3],lLf]: 

^ + ^ + 31(^-^)^ = 0. (2.9) 

Here, the asterisk denotes the complex conjugate. Using a simple point transformation 
(cf. refs. EUESD, we can convert (12. 9p into an NLS-type equation perturbed by higher 
order terms. 



2.1.3 Relation to the second flow of the matrix NLS hierarchy 



The matrix DNLS hierarchies, including the matrix Chen-Lee-Liu hierarchy, can be con- 
sidered embedded in a generalization of the matrix NLS hierarchy [23J. In the case of 
reduction (12.61) . the particular relation qr = simplifies the embedding formulas to a 
considerable extent. As a result, (12. 7p can be derived from the second nontrivial flow of 
the matrix NLS hierarchy, i.e., a matrix analogue of the (non- reduced) complex mKdV 
equation [361137] 

Qt + Qxxx ~ 3Q X RQ — 3QRQ X = O, 
Rt + Rxxx — 3R X QR — 3RQR X = O, 

through a simple, but not ultralocal, reduction. Indeed, if we set 



(2.10) 



Q 



(Ml, . . . ,U M ) 



R = CQ 



T 



V Yl!k=l ^MkUk,x J 



(2.11) 
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in terms of an antisymmetric constant matrix C, we obtain the following relations: 

Q X R = 0, QR = — Cj k (uj^ x u k — uju k)X ), QR X = (QR) X . 

l<j<k<M 



Thus, the two matrix equations f 1 2 . 1 1) simply collapse to form a single vector equation 

There exist a few advantages in regarding (12. 7p as a reduced form of (I2.10p . First, an 
infinite set of conservation laws for (I2.10p can be constructed systematically and rather 
easily using a recursive formula based on the Lax pair [3B] ■ Thus, we can obtain the con- 
servation laws for (12. 7p from those for (I2.10p through the reduction (12. lip . The conserved 
densities of the first two ranks obtained in this manner are given by 



Uj :X Uk (j 7^ ^Oj ^ C jk{Uj^ xx U k ^ x Uj^Uk,; 

l<j<k<M 



^ ] Cjk(uj :X Uk — UjUk y 



l<j<k<M 



which are derived from the conserved densities RQ and tr (Q X R X + QRQR) of (I2.1(jp . 
respectively. Second, considering a similar reduction for a space-discrete analogue of 
(I2.10p [23] . we can obtain an integrable semi- discretization of the continuous system (12. 7p . 
We employ this approach in subsection 13. II since none of the (semi-) discrete analogues of 
the third-order matrix Chen-Lee-Liu system (12. 2p . which are integrable and permit a 
reduction like (I2.6p . are known. 



2.1.4 Solution formulas 

An approach based on the inverse scattering method [26] enables the derivation of the 
solutions of the matrix Chen-Lee-Liu hierarchy, that is, (12. ip and (12. 2 p as well as the 
higher flows, through a set of formulas (cf. ref. [39] for the case of scalar variables) 



r 



K(x, x) 
K(x, x) 



K(x,y) = F(y) + - 



K(x,y) = F{y) 



dsi 



dsi 



ds 2 K{x 1 Si)F(si + s 2 — x 
dK(x, si) - 



dF{s 2 + y 
ds 2 



ds 



dsi 



F(s 1 + s 2 - x)F(s 2 + y — x) 



(2.12a) 
(2.12b) 

V > x, 
(2.12c) 

y > x. 
(2.12d) 



Here and hereafter, the bar does not denote the complex/Hermitian conjugate in general. 
The time dependence of the functions is suppressed in formulas (12.121) . The functions F(x) 
and F(x) satisfy the corresponding linear uncoupled system of matrix PDEs, e.g., 

OF d 2 F OF d 2 F . . 

'^ + a? = °- 'W 2 -l^ = (2 ' 13) 

for the second-order matrix Chen-Lee-Liu system (12. ip and 



dF d 3 F 
dt 3 dx 3 



O. 



dF d 3 F 

dt 3 dx 3 



O 



(2.14) 
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for the third-order matrix Chen-Lee-Liu flow (12.21) . and decay rapidly as x — > +00. Note 
that formulas (12.121) involve only equal-time quantities and depend on the time variables 
through the time evolution of F and F. 

The reduction (12.61) is achieved at the level of the solution formulas by setting 

F(x, t) = (fx, f 2 , . . . , f M ){x, t) =: f{x, t), F{x, t) = 2iCF(x, t f . 

With this reduction, the set of formulas (12.121) is reduced to a compact form. Thus, the 
solutions to the coupled derivative mKdV equations (12.71) . decaying as x — > +00, can be 
constructed from those of the linear vector PDE f t + f xxx = through the formula 

u(x,i) = k(x,x;t), (2.15a) 

uf \ *f \ r a r a uf , T df(s 2 + y-x) 

k(x,y) = f(y) - J dsi / ds 2 k(x,si)Cf{s 1 + s 2 -x) — , y > x. 

(2.15b) 

Here, u = (u\,u 2 , • • • ,«m) and k(x,y) are M-component row vectors. Note that for a 
given /, the integral equation (12.15bj) . referred to as the Gel'fand-Levitan-Marchenko 
type, is linear for unknown k. We assume that both u and the eigenf unctions bound in 
the potential u (cf. (I2.5P ) should also decay as x — > —00. Then, we can verify that the 
number of distinct exponential functions comprising f(x,t) has to be even. Substituting 

f(x, t) = a ie iAl * +iA ?* + a 2 e iX2X+iX i\ Im A, > (j = 1, 2), Ai ^ A 2 , (a x C, a 2 ) := ai Ca T 2 ± 0, 
k{x,y;t) = kxix,^ 1 ^' 1 ^ + k 2 {x, t)e iX2y+iX ^ 

into (I2.15P and solving it with respect to ki and k 2 , we obtain the "unrefined" one-soliton 
solution of system (12. 7\f , 



u(x,t) = fc!(a:,t)e iAlx+iA ?' + k 2 (x,t)e lX2X+ ' lX ^ 



(2.16) 



Note that the denominator in the above expression may become zero for certain values of 
x and t. By introducing a new parametrization, 

i(Ai - A 2 ) ^ = _ e _ 25 ^ ^ c ^ a =: 2e s b ^ fl2 = . 2e - 5 b 2l 



2(Ai + A 2 ) 2 
(I2.16P can be rewritten as 

5 e i(Ai-A 2 )x+i(A?-A3)t — I — 5 -i(Ai-A 2 )x-i(Af-A3)t 

u(x,t) = — r- : -— = , (2.17) 

K 1 cosh[i(Ai + A 2 )a; + i(A? + A|)t-5] V ; 

with the condition -2i(A x - A 2 )(&iC, b 2 ) = (Ai + A 2 ) 2 . The "one-soliton" solution f[2~TTj) 
resembles the soliton solution of the vector mKdV equation (11.21) [38] . 

Cie ^(^i-^2)x+^{Xl-Xl)t _j_ c e -|(Ai-A 2 )a;-i(Af-A|)t 

q ^ = cosh [i(Ai + A 2 )a; + |(Af + X 3 2 )t - 5] ' 
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under the conditions (ci,Ci) = (c 2 ,c 2 ) = and — 8(ci,c 2 ) = (Ai + A 2 ) 2 . Moreover, if 
we impose the "reality conditions" A 2 = —A* and e 25 ^ ffi<o, (12. 17[) provides the bright 
one-soliton solution of system (12.71) that behaves regularly for real x and t. With the 
parametrization Ai = £i + ir/i, A 2 = — A^ = — £i + ir]i ^ 0, r/i > 0), it reads as 

" (x ' t) = cosh fas + ,fr(3tf-tf)t + ,S] ' ^ C ^> = ^' Im(25) ^ vr (mod 

Similarly, the general iV-soliton solution before imposing the "reality conditions" is 
obtained by substituting the expressions 

27V 



f(x, t) = J2 aie iXlX+iXft , Im A, > 0, A, ^ A fc if j ? k, 
i=i 

2N 



i=i 

into (I2.15P and solving the resulting linear algebraic system for ki (I = 1,2,..., 2N) using 
a matrix inversion. In addition, conditions such as 

Pfaffian of ( ^ 

have to be imposed for the solution to decay as x — > — oo and to behave properly as 
solitons. The details will be published elsewhere. It should be noted that the soliton 
solutions thus obtained involve more free parameters, and hence, are more general than 
the previously known solutions [T2|[T3]. 

2.2 Coupled derivative mKdV equations of Kaup— Newell type 

2.2.1 Lax pair for the second flow of the matrix Kaup— Newell hierarchy 

Another DNLS equation iq t2 + q xx ± i(\q\ 2 q) x = 0, which is often referred to as the Kaup- 
Newell equation [ID], permits an integrable matrix generalization [T6 | [T7 ][30lHT] 

i% + qxx ~ i(qrq) x = O, 
in 2 - r xx - i(rqr) x = 0. 

Here, q and r are l\ x l 2 and Z 2 x li matrices, respectively. In this paper, we are more con- 
cerned with the next higher flow in the matrix Kaup-Newell hierarchy that commutes with 
the first nontrivial flow H 2 . 1 8 1) . It is written as (up to a scaling of the time variable) [16] 

qt 3 + qxxx ~ i^{q x rq + qrq x ) x - ^{qrqrq) x = O, 

3 3 (2 ' 19) 

r ta + r xxx + i-(r x qr + rqr x ) x - -(rqrqr) x = O. 
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Note that the lower indices of t are omitted in the following. The Lax pair for ( 12 . 1 9 j) is 
given by 



u = K 2 



-h 



+ C 



(2.20a) 



V = iC 

+ c 3 

+ iC : 

+ c 



-Ah 



c 



4q 

4r 

2iq x + 2qrq 



4/ 2 

— 2ir x + 2rqr 
-i(q x r - qr x ) - \{qr) 2 



+ K 4 



—2qr 



2rq 



K r Qx ~ r x q) + Urq) 2 



-r xx ~ if + rqr x ) + |(rg) 2 r 



-Qxx + iUq x rq + qrq x ) + Hqr) 2 q 



(2.20b) 



Substituting the Lax pair (I2.20p in the zero-curvature condition (12. 4p . we obtain the 
third-order matrix Kaup-Newell system ( 12.191) without any contradiction or additional 
constraint. 



2.2.2 Reduction 



Both the first flow (12.181) and the second flow (12.191) of the matrix Kaup-Newell hierarchy 
permit the reduction of the Hermitian conjugation [33] r = v4 1 g^y4 2) where A\ and A 2 are 
constant Hermitian matrices: A\ = A iy A i<t = A ifX = O. Moreover, unlike the first flow 
(12.181) . the second flow (12.191) allows an interesting reduction such that r is identically 
equal to Cq T , where C is an antisymmetric constant matrix: C T = —C, Ct = C x = O. For 
the reduced system to assume a concise form without the imaginary unit and fractions, 
the following vector reduction (h = 1, h — M) is considered: 



q = (ui, . . . , u M ), r = 2\Cq T = 2i 



V Sfell ^MkUk J 



(2.21) 



which also implies the simple relation qr = 0. System (I2.19P is then reduced to another 
system of coupled derivative mKdV equations: 



dui d z Ui d 
dt dx 3 dx 



X<j<k<M 



duj_ 
dx 



U k - Uj 



duk 
dx 



1,2,. ..,M. 



(2.22) 



To the best of the author's knowledge, this system was reported for the first time in [15] . In 
contrast to (12.71) . the partial differentiation with respect to x acts on the entire nonlinear 
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term in (|2.22p . The conserved densities of the first three ranks are given by 
Ui (i = l,2,...,M), 



l<j<k<M 

^ C 'jk(.'U>j i xx'U>k,x Uj^ x Uk t xx^) 2 
l<j<k<M 



/ J C jk (u jtX U k - UjU k ,x) 
X<j<k<M 



Using the vector notation, ( 12. 22ft can also be written as 

u t + u xxx + 3((u x C, u)u) x = 0, C T = -C. 



The Lax pair for this reduced system is obtained by substituting ( I2.2ip into q and r in 
( I2.20P ; using a simple gauge transformation, it can be rewritten in the form 



(2.23a) 



u = 


—A u 
\Cu T 


1 






r a 3 + 


2X(u x C, u) 


—X 2 u + \u x — u xx — 3(u x C, u) u 


V = 


-X 3 Cu T - 


- A 2 C< 

- 3\(u x C, u)Cu T 


\ 2 Cu T u + XC(u x u — u T u x ) 



(2.23b) 

where u is a row vector and A := 2i£ 2 . In the simplest nontrivial case of M = 2, setting 

Cxi = i, mi = i/j, u 2 = if)*, 
the following single equation is obtained: 

A + ^xxx + 3i U^ x r - MM r = 0- ( 2 - 24 ) 



Using a simple point transformation (cf. refs. EHESj), (I2.24p can be converted into an 
NLS-type equation perturbed by higher order terms. 

2.2.3 Solution formulas 

A set of formulas for the solutions of the matrix Kaup-Newell hierarchy, decaying as 
x — > +oo, can be derived by exploiting its relationship [16j[T7] with the matrix Chen- 
Lee-Liu hierarchy studied in subsection 12.11 It is written in the form of the product of 
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two matrices as follows [26J (cf. ref. |39]for the case of scalar variables): 

q = K(x,x)[I + L(x 1 x)], (2.25a) 

r=[I + L(x, x)Y l K{x, x), (2.25b) 

K(x,y) = F(y)+ 1 - J°° d Sl J°° ds 2 K(x, Sl )F( 8l + s 2 - x) ~ %) , V>x, 

(2.25c) 

• roo poo dK(x Si) - 

K(x, y) = F{y) --J d Sl J ds 2 ^ 1J F( Sl + s 2 - x)F{s 2 +y-x), y>x, 



i f°° 

L ( x ,y) = -2j dsF(s)F(s + y-x) 



(2.25d) 



j roo poo OL(x s ) 
+ - / dsi ds 2 K ' 1J F( Sl + s 2 - x)F(s 2 +y-x), y>x. (2.25e) 
2 J x J x osx 

The time dependence of the functions as well as the index of the unit matrix I to indicate 
its size is suppressed in the above formulas. Here, F(x) and F(x) satisfy the corresponding 
linear uncoupled system of matrix PDEs, e.g., (I2.13P for the matrix Kaup-Newell system 
(12.181) and (12.141) for the third-order flow (I2.19p . and decay rapidly as x — > +oo. This set 
of formulas is indeed useful and sufficient for constructing explicit solutions of the matrix 
Kaup-Newell hierarchy. However, its shortcoming is that one of the most important 
properties of the hierarchy has not been incorporated, that is, the hierarchy allows the 
introduction of the potential variables q =: q x and r =: r x , and can be reformulated in 
terms of q and f. Consequently, the elementary function solutions of the matrix Kaup- 
Newell hierarchy can be expressed as the partial x-derivatives of elementary functions; 
this fact has passed unnoticed in the existing literature (see refs. [16l[22] for some indirect 
results). To bridge this gap, we propose a new set of solution formulas that accurately 
reflects the feasibility of potentiation for the matrix Kaup-Newell hierarchy; it assumes 
the following form (^-derivatives of single quantities [TBII2^] ): 

dhC(x,x) 

q = -^-L, (2.26a) 

r = 8 ^ X \ (2.26b) 

r°° i r°° r°° dlC(x si) 

JC(x,y) = -l ds 1 F(s 1 ) + -J dsxj ds 2 ^ 1J F( Sl + s 2 - x)F{s 2 + y - x) 

% i f°° j r A dlC(x,s 1 )dG(s 1 + s 2 -x)dG(s 2 + y-x) 
-G(y) + -l ds,^ ds 2 — , y>x, 

(2.26c) 

iC(x, V) = - f d*i F( Sl ) d Sl 1°° ds 2 9 ^ d X J^ F( Sl + s 2 - x)F{s 2 + y-x) 

s i f°° a Ha 9iC(x,s 1 )dG(s 1 + s 2 -x)dG(s 2 +y-x) 

= G{y) -2L dsl L dS2 ^ dy ' y ~ x - 

(2.26d) 
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Here, the matrices G and G are the primitive functions of F and F, respectively, that also 
decay as x — > +00, that is, G(x) := — J°° F(y)dy and G(x) := — j°° F(y)dy. Note that 
d/dx in f !2.26ap and ( 12. 2 6b j) denotes the partial differentiation with respect to x, while 
all the time variables are fixed. The two sets of formulas ( I2.25P and (I2.26P are indeed 
equivalent, though it is a highly nontrivial task to verify it directly. 

Similarly to the case of the matrix Chen-Lee-Liu hierarchy, the reduction (I2.2ip can 
be realized by setting 



G(x, t) = (g 1 , g 2 ,..., 9m){x, t) =: g(x, t), G(x, t) = 2iCG(x, t) T . 



(2.27) 



Thus, the set of formulas (I2.26P is reduced to a compact form. In particular, the so- 
lutions to the coupled derivative mKdV equations ( 12.22(1 . decaying as x — > +00, can be 
constructed from those of the linear vector PDE g t + g xxx = through the formula 



d 

u(x,t) = —k(x,x;t), 

poo 

Hx,y) = g(y) - / dsi 



00 ^ dk(x, s^ ^dgjsi + s 2 - x) T dg(s 2 +y-x) 



dsi 



dy 



(2.28a) 

y > x. 
(2.28b) 



Here, u = (ui,u 2 , ■ ■ ■ ,Um) an d k(x,y) are M-component row vectors. Substituting the 
expressions 



g(x,t) = aie iAia;+iA ^ + a 2 e iA2a;+iAi * 



^, ImXj > (j = 1,2), Ai ^ A 2 , (a x C, o 2 ) := aiCaJ ^ 0, 
k{x, y; t) = k x (x, t) e iA ^ +iA ?' + k 2 (x, t) e iA2y+iA '* 



into (I2.28P and solving it with respect to ki and k 2 , we obtain the "unrefined" one-soliton 
solution of system (12.221) . 



u(x, t) 



d_ 

dx 
d_ 

dx 



ki(x, t)e 



iAix+iAf t 



+ k 2 (x,t)e 



die 



+ a 2 e 



iA2a;+iAot 



1 + 



iAiA2(Ai— A2) 
2(A!+A 2 ) 2 



OjC, a 2 )e i ( Al+A2 )^+ i ( A f+ A i) i 



(2.29) 



Note that the above solution also decays as x — > —00, but may have singularities for 
certain values of x and t. By introducing a new parametrization, 



[XlX2{Xl ~ X2 ^ ai C,a 2 ) =:e~ 2S 



2(Ax + A 2 ) 2 
(12 . 29(1 can be rewritten as 



(5 e C), ai 



--: 2e~ s b ll a 2 =: 2e- 5 b 2 , 



u(x, t) 



d_ 

dx 



5 ie |(Ai-A 2 )a : +i(Af-A|)t + 5 2e -^(Ai-A 2 )a;-|(Af-Al)t 

cosh [i(Ai + A 2 )x + i(A? + \\)t - S] 



(2.30) 



with the condition 2iAiA 2 (Ai — X 2 )(biC,b 2 ) = (Ai + A 2 ) 2 . Moreover, if we impose the 
"reality conditions" A 2 = — A^ and e 2S M<o, ( 12.30(1 provides the bright one-soliton solu- 
tion of system ( 12.221) that behaves regularly for real x and t. With the parametrization 



15 



Ai = £1 + ir/i, A 2 = —XI = — £1 + ir)i (£1 ^ 0, 771 > 0), it reads as 



u(x, £) 



9 J b ie i«iz+i€i(d 2 -3r? 2 )i + 5 2e -i6^-i6(?i 2 -3r? 2 )t 



<9x I cosh [rjix + ?7i(3^ — r]f)t + 5] 




bi cosh [ijxx + 771(3^1 - r/i)t + 5 + i<p] e 



,i?lz+i«i(£i-3»7 2 )* 



b 2 cosh [7/1 a; + 771 (3£i - r/j^t + 5 - i<p] e 



^l^-i«i(?i 2 -3r?' 




with the conditions i£i(biC, 6 2 ) = + »7i)> exp(iy2) := (£1 + ir)i)/y/g[ + rjf, and 

Im(25) ^ 7T (mod27r). 

2.3 Massive Thirring-like model with symplectic invariance 
2.3.1 Derivation 

Using our approach, it is possible to obtain not only evolutionary systems but also non- 
evolutionary systems with symplectic invariance. To demonstrate this, let us consider the 
application of the same type of reduction as described in subsections 12.11 and 12.21 to the 
first negative flows of the matrix DNLS hierarchies. The first negative flow of the matrix 
Chen-Lee-Liu hierarchy reads as [T6tlTT] 



while the first negative flow of the matrix Kaup-Newell hierarchy reads as [T6|IPT] 



These are the matrix generalizations of massive Thirring-type models [T8H22] . where m 
denotes an arbitrary nonzero constant responsible for the mass terms. The massless limit 
m — > is not considered in this paper; hereafter, the value of m is set as 1, without loss 
of generality. This is easily achieved by rescaling 0, x, and d T . The Lax pairs for systems 
(12.311) and (12.321) are already known [T6UTT] , and thus, they have been omitted here. Both 
the systems permit the reduction r = Cq T , \ = —C<f) T , C T = —C, which is the natural 
extension of the reduction r = Cq T , C T = —C for the positive flows. 

Let us consider a further reduction to the case of row/column vector variables. For 
the reduced systems to assume a concise form without the imaginary unit and fractions, 



iq T + mcj) - §0X<7 
ir T — mx + \t4>X 
i(f) x + mq — ^<f>rq 
iXx -mr + \rqx 



O, 

o. 

o, 
o, 



(2.31) 



iq T +m(j)- \(qx<P + <Pxq) = O, 
ir T -mx + |(r0x + X<pr) = O, 
i(j) x + mq = O, 
iXx — Tnr = O. 



(2.32) 
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we consider the following vector reduction: 

q = Oi, . . . , m m ), r = 2iCq T = 2i 



( Y,k=iCikV k \ 



i(t>i, . . . , v M ), X = -2iC0 T = 2 



Then, because (fix — 0, system (I2.3ip reduces to 



V Y.iLiCukVk J 



d 2 Uj 
dxdr 



+ Ui- 



E c 



_l<j<k<M 

system (I2.32p collapses to 

d 2 v,; 



drdx 



^2 C i k 

l<j<k<M 



dr 



Uk — u 



du k \ 
dr J 



Ul = 0, i = l,2,...,Af; 



(2.33) 



V dx 3 dx 



Vi = 0, i = l,2 



,...,M, 



(2.34) 



where Ui = dvi/dx, % = 1, 2, . . . , M. These are considered non-evolutionary symmetries 
of system (12.71) and system (12.221) . respectively. It can be easily observed that (12.331) 
and (12.341) are equivalent up to a simple interchange of the variables. Nonetheless, the 
commutativity of the negative and positive flows in each hierarchy itself is of significance. 
The Lax pair for (1231) is given by (cf. (jZgSJ) ) 



U 



—X v x 
XCvl O 



V 



A A 

Cv T -Cv T v 



(2.35) 



where v = (vi,v 2 , ■ ■ ■ ,%). In addition, it may be noted that (I2.34|) can be rewritten as 
a system for Ui{= Vi )X ). Indeed, ( 12.341) implies the relation 



\^ n dui 

l<i,l<M 

from which we obtain 



E c 



jk 



l<j<k<M 



/dVj_ 

V dx 



V k - Vj 



dv k \ 
dx J 



Cuviui = o, 



Ki,KM 



E fdvj dv k 



l<j<k<M 



dx 



(2.36) 



under appropriate boundary conditions. Substituting (12.361) in (12.341) with a simple oper- 
ation, we arrive at a closed PDE system for Uj, 



d_ 

dx 



dui 
> i 

'dr 



+ Mi = 0, i = l,2,...,M. (2.37) 
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2.3.2 Solution formulas 

A set of formulas for the solutions of the matrix massive Thirring model ( 12 . 3 1 f) with 
m = 1, decaying asx-)> +oo, is given by (12.121) supplemented with the following: 

-i</>= J(x,x), (2.38a) 
i X = J(x,x), (2.38b) 

J(x,y) = -J dsF(s)-^J dsi J ds 2 J(x, s ~~~^~^ — —F(s 2 + y - x), 

y>x, (2.38c) 

1( \ f°°A w< \ 1 f°°A f°°A dJ(x, Sl )dF(s 1 + s 2 -x) [°° 

J(x,y) = - dsF(s)-- ds 1 ds 2 — ^ ~ / ds 3 F(s 3 ), 

Jy 1 Jx Jx ° S 1 ° S 2 J S2+y - x 

y>x. (2.38d) 

Similarly, a set of formulas for the solutions of the matrix Kaup-Newell-type Thirring 
model (I2.32p with m = 1, decaying as x — > +oo, is given by (I2.26P together with the 
following two relations: 

-i<f> = K.(x,x), (2.39a) 
iX = 1C{x,x). (2.39b) 

In both cases, the linear matrix PDEs satisfied by F and F are given by 

d 2 F - d 2 F 

+ F = 0, —— + F = 0, 



drdx drdx 

and the same relation applies for their primitive functions G and G. Applying the same 
reduction as that in the positive flow case (cf. (I2.27P ) to the formulas (I2.26P and (I2.39p . 
we obtain the solution formula for system ( 12.341) . 

v(x, t) = k{x, x;t), (2.40a) 

u x ( , r A [°° A dk ^ g i) ^ d 9(si + s 2 - x) T dg(s 2 + y-x) 

k{x,y)=g(y)- ds 1 ds 2 — C , y>x. 

Jx Jx os i os 2 Oy 

(2.40b) 

Here, v = (vi, v 2 , . . . , %), and g solves the linear vector PDE g XT + g = 0. Using the 
formula (I2.40p . we can construct the soliton solutions of system (12.341) in a manner similar 
to that in the positive flow case. In particular, the one-soliton solution of ( I2.34[) is given 
by 

v(x,t) = . , i^i(6iC,62)- 



cosh h] lX - ^jjT + 5) ?i + Vi 

The additional conditions (i 6l - {0}, 771 > 0, and e 2<5 ^ IR<o guarantee that this solution 
is regular for real x and r. 
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3 Integrable discretizations 



In this section, we propose space discretizations of systems (12.71) . (12. 22ft . and (12. 34ft while 
retaining both the integrability and symplectic invariance in the continuous case. Using a 
tricky transformation peculiar to the discrete case, we also obtain a novel integrable semi- 
discretization of the Manakov model (jl.ip . One-soliton solutions are obtained by solving 
discrete integral (linear summation) equations of the Gel'fand-Levitan-Marchenko type. 



3.1 System of coupled derivative mKdV equations (12.71) 

A space-discrete analogue of the (non-reduced) matrix complex mKdV equation (I2.10p 
that is "maximally symmetric" with respect to reduction of the matrix variables is given 
by [23] 

Qn,t + (h — Qn+lRn+l) l Qn+l + (h _ QnRn) 1 Qn 

(h — Q n R n+ i) l Q n — (h — Q n -\R n ) Q n -i = O, 

(3.1) 



Rn,t + Rn+l (h — QnRn+1 



Rn{I\ — Qn-lRn 



,-1 



o. 



— Rn(h — QnRn) 1 — Rn-l(h — Qn-lRn-l) 1 

Here, "maximally symmetric" implies that ( 13. ip possibly permits all the interesting reduc- 
tions corresponding to those in the continuous case. System (13.11) possesses the following 
Lax pair: 



Zh + QnRn \Qn 

-(Z + I) Rn \h 







Qn 









J 2 . 





h o 

{z + \)R n \h 



-( z * + js+2)h 

+ (l + £) (h-Qn-lRn)- 1 


— {h — QnRn) 
— -^l{h — Qn- 


Qn 

-\Rn)~ 1 Qn~l 


{z 2 + \){h-R n . 1 Q ri . l )- l R n ^ l 

+ (1 + jl) {h ~ RnQn-l)- l Rn 


-(i + i) (h 


RnQn—l) 



(3.2a) 



(3.2b) 



Here, z is the spectral parameter independent of n and t; h and h are the l\ x / x and 
/ 2 x l 2 unit matrices, respectively. Indeed, substituting (13. 2p in (a space-discrete version 
of) the zero-curvature condition [12] 

L n , t + L n M n - M n+1 L n = O, (3.3) 

which is the compatibility condition for the overdetermined linear equations \l/ n +i = L n ^ n 
and ty ni t = M n ^ n , we obtain the space-discrete system (13.11) . To achieve a discrete coun- 
terpart of the reduction (12.111) . we first assume the following relations between Q n (n G Z) 
and R n (n G Z): 

R n = P n -P„„ 1 , Q m Pn + Q n P m = 0, Vm,neZ. (3.4) 
If these are satisfied, then system (13 .ip is reduced to the form 

Qn,t + (h + Qn+lP-n) (Qn+1 ~ Qn) + (h + QnPn-l) 1 (Qn ~ Qn-l) = O, 

[Pn,t + {Pn+1 - Pn){h ~ QnPn+l)^ + {P n ~ P n -l){h ~ Qn-lPn)^] ~ [n -> 71 - l] = O. 

(3.5) 
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The following choice of Q n and P n automatically satisfies (J3T4J) : 

Qn \ • • ■ i ^ ) ) -f n CQ n 



c 



r 



-C. 



(3.6) 



Substituting (13. 6p in (13. 5p . we obtain an integrable semi- discretization of the coupled 
derivative mKdV equations (12. 7p . 



l<j<k<M 



+ 



(i) _ (»') 
Un U n-1 



U n' U n+l) 1+ E/ — U n~l U n^) 

l<j<k<M 

Z = 1,2,...,M. (3.7) 

The Lax pair for the reduced system (13.71) is obtained by substituting (I3.6P into Q n and 
R n {= P n — Pn-i) in (13.21) . Moreover, using a gauge transformation, we can restore the 
ultralocality of the spatial Lax matrix L n ; thus, the Lax pair can be written in the form 

/i (fi + l)u n 

(ji-l)Cul / + (// + !) Cu^u n 



0. 



(3.8a) 



-1 



1 + (u n C, u n -i) 







_ 1 


(/i+i 


) l*n + (l + W n _i 


(A* 
+ 


"I) 

o- 


- i) c < 


(ah 

+( 


- 1) CwJ_ x w n 

1 + C«n«n-1 



(3.8b) 



where u n = (m^ ; , u^, • • • , u[ M ') is a row vector and fi : = z 2 is the "new" spectral parame- 
ter. It is easy to see that the discrete eigenvalue problem *f? n +i = L n ^/ n with (I3.8al) reduces 
to the continuous eigenvalue problem ty x = with (I2.8al) in a suitable continuous limit 
(cf. ref.HiD. 



3.2 Manakov model (ECO 



Let us consider the canonical case of the coupling constants (cf. the introduction) in (13. 7p . 
namely, C 2 j-i2k = -C 2k 2j-i = S jk , C 2 j-i2k-i = C 2 j2k = 0, and M = 2m, and change the 
dependent variables as follows: 

u (2j-d = . ( _i)v*gW u m =: in+V^rW, j = 1, 2, 

Then, system (13. 7p is converted into the following form: 



m. 



„(i) :„(i) 



+ 



- 2gW = 0, 



fc=i 



fe=i 



. <9r 



(?) 



' n+1 ^ 



r (j) _:„(*) 
'n-1 



<9f 



+ 2ri j) = 0, 

1 + E feffirf' + rf^fflO 1 + E GMi + ) 

J = l,2,...,m. 



fc=i 



fe=i 
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By further imposing the reduction r$ = ■> °j — =tl> we obtain a new integrable 

semi-discretization of the system of coupled NLS equations. The simplest case Tn = —qtP* 
(j = 1, 2, . . . , m) provides a space-discrete analogue of the £7(m)-invariant Manakov model 

(HH), 

i Qn+1 ~ 1( ?n | Qn-1 + l Qn r, _ q /o q\ 

1<ln " 1 - <<Zn + l, Q*n) - (Qn, Qn + l) 1 " (Qn, Q*n-l) ~ tin-l, Qn) ^ ~ ' 1 ' J 

where q n = (q£\ q n 2 \ . . . , gi m) ) . 



3.3 Solutions to system ( 13.71 ) 



A set of formulas for the solutions of the space-discrete matrix complex mKdV system 
(13. ip as well as its commuting flows, which tend to zero as n — > +00, is given by 

Q n = K(n,n) + K(n + l,n+ 1), (3.10a) 

R n = K{n,n) + K{n + l,n + l), (3.10b) 

00 00 

K(n, m) = G(m) + \- K ( n > n + •?') + K ( n ' n + 3 + l )\ l G ( n + j + k + 1) + G(n + j + k + 2)] 

x [G(m + k) + G(m + fc + 1)] , m > n, (3.10c) 



00 00 



K(n, m) = -G(m) + ^ [ K ( n > n + o) + n + J + 1)] + 3 + k ) + G ( n +j + k + l 

x [G(m + k) + G(m + k + 1)] , m > n. (3.10d) 

Here, the functions G(n) and G(n) satisfy the corresponding linear uncoupled system of 
matrix differential-difference equations, e.g., 

+ G(n + 1) - G(n - 1) = O, + C(n + 1) - G(n - 1) = O 

for the flow ( 13. ip . and decay rapidly as n — >■ +00. The reduction given by ( 13. 4ft and ( 13. 6 p 
is achieved at the level of the solution formulas by setting 

G{n) = ( 9l ,g 2 , . . .,g M )(n) =: g(n), G(n) = -C[G(n) - G(n - l)f. 

With this reduction, the set of solution formulas (I3.10p is reduced to a compact form. 
Thus, the solutions to the semi-discrete coupled derivative mKdV equations (13. 7p . decay- 
ing asn-f +00, can be constructed from those of the linear vector differential-difference 
equation dg{n)/dt + g(n + 1) — g(n — 1) = through the formula 

u n (t) = k(n, n; t) + k(n + 1, n + 1; t), (3.11a) 



00 00 



k(n, m) = g(m) + ^2^2 [ k ( n > n + + fc (™, n + 3 + ^)] G [d( n + 3 + ~~ fl 1 ^ + J + 1 + 2 )] 

i=o 1=0 

x [g(m + l)+g(m + l + l)], m>n. (3.11b) 
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Here, u n = (v,n , Un \ ■ ■ ■ , Un) and k(n, m) are M-component row vectors. Substituting 
the expressions 

g(n,t) = a^e^-^ 1 + a 2 ^ n e^-^\ \ H \ > 1 (j = 1,2), ^ ^ /i 2 , ( ai C,a 2 ) 
k(n,m;t) = fci(n, t)^ m e^'^ 1]t + fc 2 (n, t)^™^-^ 

into (13. lip and solving it with respect to k\ and k 2 , we obtain the "unrefined" one-soliton 
solution of system (13. 7h in the additive form 

«n(f)=P„(*)+P„+i(*), (3-12) 

where p n (t) is given by 

p n (f) = fe^n,*)/^*" 1 -"^* + fc 2 (n,t)^J n e (w -'^ 1) * 
oi/ipe^ 1- ^ 1 ) 1 + a 2 fi 2 n e < - fl2 -^ 1)t 

(3.13) 



In fact, the row vector p n (t) itself satisfies the nonlinear differential-difference equation 
(cf. (EID and (EH), 

^Pn _| Pn+l ~ Pn-l _ q /g ^\ 

^ 1 + (Pn+lC> Pn> + (PnC; Pn-l) + (p n +l C > Pn-l) 

Note that the denominator in the expression (13.131) may become zero for certain values of 
n and t. By introducing a new parametrization, 

(pi ~/i 2 )(l + /ii)(l +/i 2 ) 25 

(aiC,a 2 )=:e (o G C), 

(1 - /ii/i 2 ) 2 

ftl =: 2e- 5 6 1; a 2 =: 2e~ 5 b 2 , ^ = e a ~ ip , /i 2 = e a+i/3 , 



(I3.13P can be rewritten as 

^^gi/3n— 2i(cosh a sin/3)t _|_ ^^g— i/3n+2i(cosh a sin /3)t 

Pntt) = — TT of ■ C T xi ' ( 3 - 15 ) 

cosh [an — 2(smh a cos p)t + oj 

with the condition — 4i(cosh a + cos (3) sin /3 (&iC, 6 2 ) = (sinha) 2 . The "one-soliton" so- 
lution (I3.15P of (I3.14p resembles the soliton solution of the semi-discrete vector mKdV 
equation (or, the vector modified Volterra lattice) dq n /dt = (l + (q n , q n ))(q n+ i — q n -i) 
given by 

_ i/3n+2i(cosh a sin/3)t _i_ _ — i/3n— 2i(cosh a sin f3)t 

q (t) ^2e 
n cosh [an + 2(sinh a cos (3)t + 5} 

under the conditions (c\,Ci) = (c 2 , c 2 ) = and 2(ci,c 2 ) = (sinha) 2 . Moreover, if the 
"reality conditions" a > 0, < (3 < n (or — n < (3 < 0), and e 2<5 ^ R<o are imposed, (I3.15P 
provides the bright one-soliton solution that is indeed regular for real n and t. Owing 
to the discrete nature of the space variable n, there can exist other cases wherein the 
solution is regular, e.g., a > 0, (3 = ir/2, and e 2an + 25 — 1 ; \/n. 
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3.4 System of coupled derivative mKdV equations ( 12.22 



In this subsection and the subsequent subsections, we use the forward difference operator 
A n to indicate 

^nfn+j - = fn+j+1 fn+j- 

It can be recalled that an integrable semi-discretization of the third-order matrix Kaup- 
Newell system (12.191) is given by [23] 



q n>t + A n [(h - q n r n ) l q n + (h + q n ~ir n ) 1 q n -i] = O, 
r njt + A n [r n {I x + q n -ir n )~ l + r„_i(/i - <7„_ir n _i) _1 ] = O. 



This system possesses the following Lax pair: 



zh ~(z-\) q n 



M n 



- (1 - jf) (h + qn-xTnY 1 







q n 









J 2 . 





h 



(3.16) 



(z 2 - 1) (I 2 - r n _ig n _i) Vn-i 
+ (l-pr) ( / 2 + r„g„-i)- 1 r n 



[h - q n r n ) 1 q n 

\{h + <7 n -i7" n )~V-i 



(l - 4,) (J 2 + r n g n _ 



(3.17a) 



(3.17b) 



Indeed, the substitution of (13. 17ft in the zero-curvature condition (13. 3ft gives the space- 
discrete system (13 . 16f) . As system (12. 19ft permits the reduction r oc Cq T , C T = —C, so 
system (I3.16P allows the corresponding reduction r n = Cq T _ L , C T = —C. In particular, 

n 2 

considering the vector reduction 



q n = (u 



(i) v (Mh 

n > • • • ) w n ) 



cq: 



c 1 



-c, 



we obtain an integrable semi-discretization of the coupled derivative mKdV equations 

dm, 



du 



(0 



Of 



+ A r 



(0 



i<i<fc<M 



+ 



■it 





n-l 



l<j<k<M 



1,2,..., M. (3.18) 



This space difference scheme depends on the five points: n, n ± ~, and nil. In fact, 
we can derive a simpler, three-point difference scheme for (I2.22p from the same matrix 
system (13.161) . For this purpose, we consider the vector reduction 



C(g n + g„_i) T , C T = -a 



(3.19) 
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Then, considering the relation qnCq^ + qmCq^ = 0, we can observe that (13.161) is reduced 
to an alternative lattice version of (|2.22p . 



du 



Of 



(i) , (*) 



i<j<fc<M 



0. 



1,2, 



M. 



(3.20) 



System (I3.20p resembles the integrable semi-discretization of the vector third-order Heisen- 
berg ferromagnet model (jl.3p (see (2.22) in ref. |T| or (5.25) in ref. HI]). The Lax pair for 
f)3.20p is obtained by substituting (I3.19P into q n and r n in (I3.17p . Moreover, using a gauge 
transformation, we can restore the ultralocality of the L n -matrix; thus, the Lax pair can 
be written in the form 



-H + l)Ct£ I 



(ji + 1) u n 

{jM - 1) Cu T n U n 



M„ 



1 - (u n C, n n _i) 





- (// + 1) w n - 




«n-l 


Oi - 1) Ctd 
+ ( x - c < 


(fi - 1) Cu^Kn 

- (l - ±) c<«„ 


-1 



(3.21a) 



(3.21b) 



where u. 



(1) 



and // := z . The discrete eigenvalue problem ^ n+ i = L n ^ n 
with ( 13. 2 lap reduces to the continuous eigenvalue problem ^f x = Uty with (|2.23aj) in a 
suitable continuous limit (cf. ref. 03]). 

It should be noted that the two systems (13.71) and (13.201) are connected through the 
change of variables 

ug> (-l)»«g>, t^-t. 

This correspondence is not surprising as their respective "ancestor" systems ( 13. ip and 
( I3.16P are connected through the same type of transformation [23] . 



3.5 Massive Thirring-like model (12.341 ) 

A space discretization of the matrix massive Thirring-type model ( I2.32p . 

q n ,r + i(0n + 0n+l) + 2(g„Xn+10n + 0n+lXn+l<?n) = O, 

r n ,r ~ KXn + Xn+i) ~ 2(r„0„x„ + Xn+10n?"n) = O, 
4> n - <f) n+ i + iq n = O, 
Xn - Xn+i ~ ir n = O, 



(3.22) 



together with its Lax pair, was proposed in ref. [23 This is considered the first negative 
flow of the semi-discrete Kaup-Newell hierarchy, which contains (I3.16P as a positive flow. 



It should be noted that system f!3.22j) permits the reduction r n = Cq 



Xh 



-c<f> r 



1 ■ 

2 



C = —C. In particular, the vector reduction 



(u 
iff. 



(i) 

n ' ' 

(1) 



Cq J 



Xr. 



c 1 



-c 
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simplifies fj3.22j) to a single vector equation, 







_< , 7j( i )_„( i M + „(*) + „« _? V r-,(v {j) v {k) -v ij) v (k) ) 



-l<j<k<M 



.0 

J n+1 



- 2 



L l<j<fc<M 



= 0, 



z = 1,2,...,M, 



(3.23) 



with «n = A n f„ . This is a three-point (n, n.+ |, n+1) difference scheme for the massive 
Thirring-like model (I2.34p . and is considered a non-evolutionary potential symmetry of 
system (I3.18p . In the same manner as described in subsection 13.41 we can also obtain a 
simpler, two-point difference scheme for (I2.34p from the matrix system (13.221) . Indeed, if 
we set 



In = {u ( n\ 



u (M) ) 

i a n )i 



r„ = C(q n + g n _i) 



ID, 



(i) 



n ■)•••■> 



T 



c T = -a 



(3.24) 



and utilize the relation n C0^ + m C0^ = 0, (13.221) is reduced to an alternative space 
discretization of (I2.34p . 



— (v^ -v (i) 



l)+ V n + V n+l / 2^ °M U n+l U n V n V n+l) 



-l<j<k<M 



(4° + «2-i) = o. 



z = l,2,. ..,M, 



(3.25) 



where wi*' 1 = A n t>i^(= — v^)- This is considered a non-evolutionary symmetry of 
(the potential form of) system (13.201) . The Lax pair for (I3.25P is given by (cf. (13.211) and 

flZSSD) 



Mr, 



1 



fJL (/j, + l)(v n+1 -v n ) 

-fi + l)C (vl +l -vl) I - (fi - 1) C (vl +l - vl) (v n+1 - v n ) 



/i+1 -2(/i + l)v n 

-2^-l)Cv T n 4(v-l)Cv T n v n 



(3.26a) 
(3.26b) 



where v n = (vn\vli \ . . . , Vn )■ Along parallel lines with the continuous case (cf. subsec- 
tion [21TT]) , we can rewrite (I3.25P as a closed differential- difference system for Un ■ The 
resultant system provides an integrable semi-discretization of system (12.371) . 



3.6 Solutions to systems ( ET20T ) and ( I3T25T ) 

In analogy with the continuous case (cf. (I2.26P ). a set of formulas for the solutions of the 
space-discrete matrix Kaup-Newell system (13.161) as well as its commuting flows, which 
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tend to zero asn-> +00, can be presented in the following difference form [26J: 

q n = A n JC(n,n), (3.27a) 
r n = A n K(n,n), (3.27b) 



00 00 



K(n, m) = - F(s) + n + ^ ~ n + 3 + X )] F ( n + 3 + k + l)F(m + k) 

s=m j=0 k=0 

00 00 

= G(m) + [^( n ' n + 3) - n + J + 1)] [G(n + j + k + 1) - G(n + j + k + 2)] 

j=0 fc=0 

[G(m + k) - G(m + k + 1)] , m > n, (3.27c) 



x 



00 00 



JC(n, m) = -^y] F(s) - ^ Y 71 + ~ ^( n ' n + 3 + F ( n + J + k ) F ( m + k ) 

s=m j=0 k=0 



00 00 



= G(m) - i^ n ' n + i)~ /C ( n > n + 3 + !)] + 3 + k )~ G(n + j + k + l)] 

j=0 k=0 

x [G(m + k) - G(m + k + 1)] , m > n. (3.27d) 

Here, the functions F(n) and F(n) satisfy the corresponding linear uncoupled system of 
matrix differential-difference equations, e.g., 

?^^. + F(n + l)-F(n-l) = 0, + Fin + 1) - Fin - 1) = O (3.28) 

at at 

for the flow (13.161) . and decay rapidly asn4 +00. The matrices G(n) and G(n) are the 
"primitive functions" of F{n) and F(n), respectively, that also decay as n — > +00 and 
satisfy the same linear system, that is, G{n) := —J2^L n F ( s ) anc ^ G(n) := —J2^L n F ( s )- 
The reduction (I3.19P is realized at the level of the solution formulas ( I3.27P by setting 

G(n) = ( 9l ,g 2 , . . .,g M )(n) =: g(n), G(n) = C[G(n) + G(n - l)f. (3.29) 

In particular, the solutions to the semi-discrete coupled derivative mKdV equations ( 13.201) . 
decaying as n — > +00, can be constructed from those of the linear vector differential- 
difference equation dg(n)/dt + g(n + 1) — g(n — 1) = through the compact formula 

u n (t) = A n k(n,n;t), (3.30a) 

00 00 

k(n, m) = g(m) + Y^2 [ fc ( n ' n + ~ fe ( n > n + j + l)]C [g(n + j + I) - g(n + j + I + 2)] 

j=o 1=0 

x [g(m + I) - g[m + 1+1)], m > n. (3.30b) 

Here, u n = (un\un \ . . . , u^) and k(n, m) are M-component row vectors. Substituting 
the expressions 

g(n,t) = a^e^-^ + a 2 ^ n e^~^\ > 1 (j = 1,2), ^ ^ /i 2 , ( ai C, a 2 ) ^ 0, 

k(n, m; t) = k^n, t)^^ 1 '^ + k 2 (n, t)^™^ 2 '^ 1 
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into (I3.30P and solving it with respect to ki and fc 2 , we obtain the "unrefined" one-soliton 
solution of system f)3.20p in the difference form 

U n (t) = A n { -. ... lft .. 1 r ^ 2 j r- \ . 3.31 

Note that the denominator in the expression (13.311) may become zero for certain values of 
n and t. By introducing a new parametrization, 

(/il ~/i2)(l ~/il)(l -^ 2 ) 7 „ n n . -25/. 

(aiC ' a2) =• 6 C) ' (3.32) 

a^^e"^!, a 2 =:2e- 5 b 2 , /ii = e a " i/3 , // 2 = e a+i/3 , 

( 13.311) can be rewritten as 

^^gi/3n— 2i(cosh a sin/3)t _|_ Jj^q~ i/3«+2i(cosh a sin/3)t 



u n (t) = A n 



cosh [an — 2(sinh a cos /3)t + 5] 



^sin (Hf*) sin (^2) 

cosh [an — 2(sinh a cos /3)t + 5' + i<p] 



cosh 2 [an — 2(sinha cos/3)t + 5'] + sinh 2 (|) 

x e i/3n-2i(coshasin/3)t _ & / cogh ^ _ 2 ( s i n ha COS 0)t + 5' - l(f] e - i ^+ 2i ( cosh a sin «* | , 

(3.33) 

with the condition 4i(cosha — cos (3) sin (3 (biC, 6 2 ) = (sinha) 2 . The constant <£> on the 
right-hand side of (I3.33|) is defined as 



exp(ip) : = 



sin 



sin (te) sin (^) 



and the new shifted parameters 5', b' l7 and b' 2 are given by 5' := 5 + |r, ft'i := e^bi, and 

6 2 := e _1 2b 2 . If we impose the "reality conditions" a > 0, < f3 < n (or — 7r < /3 < 0), 
and e 2<5 ^ M <0 , (I3.33j) provides the bright one-soliton solution of fl3.20j) that is indeed 
regular for real n and t. 

A set of formulas for the solutions of the first negative flow (13. 22f) of the semi-discrete 
matrix Kaup-Newell hierarchy, decaying as n — > +oo, is completed by supplementing 
(I3.27P with the following: 

— i(p n = JC(n, n), (3.34a) 
iXn = K-(n,n). (3.34b) 

The linear uncoupled system of matrix differential-difference equations to be satisfied by 
F(n) and F(n) in this flow reads as 

dF(n) dF(n + l) — . . -, , „ dF(n) dF(n+l) . . , N _ 
+ F(n) + F(n + 1) = 0, + F(n) + F(n + 1) = O, 
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and the same relation applies for their "primitive functions" G(n) and G(n). Assuming the 
same restriction as that in the positive flow case (cf. ( 13.29)) ). we can realize the reduction 
(13.241) on the solution formulas ( I3.27P and (I3.34p . Thus, the solutions to the semi-discrete 
Thirring-like model ( I3.25p . decaying as n — > +00, can be constructed from those of the lin- 



ear vector differential-difference equation dg(n) /dr — dg(n + l)/dr + g{n) + g(n + 1) = 
through the compact formula 



v n (r) = k(n,n; r), 



(3.35a) 



00 00 



k(n, m) = g(m) + [k(n, n + j) — k(n, n + j + 1)]C [g(n + j + I) — g(n + j + I + 2 

j=o 1=0 

x[g(m + l)-g(m + l + l)], m>n. (3.35b) 



Here, v r 



Vn ,Vn 



Using the formula (I3.35p . we can construct the soli- 



ton solutions of system (I3.25P in a manner similar to that in the positive flow case. In 
particular, the "unrefined" one-soliton solution of (I3.25P is given by 



V n T 



_ Ml+l T _ M2 + 1 . 



1 - 



(M1-M2)(1-^i)(1-^2) 



(a 1 C,a 2 )/ir> 2 - n e-^ + ^ri r 



(i-/ii/i 2 ) 2 

which can be rewritten in terms of the parametrization (13.321) as 



Vnir) 



i- 



sin /3 



+ b 2 e 



-i/3n+i- 



sin f) 



cosh 



an + 



sinh a 



2sinh(^)sinh(^^) 



1^^ + ^ 



with the condition 4i(cosha — cos/3) sin/3 (b\ C, b 2 ) = (sinh a) 2 . This provides the bright 
one-soliton solution under the "reality conditions" a > 0, < (3 < n (or — n < f3 < 0), 
and e 25 ^ M<o, which is indeed regular for real n and r. 



4 Concluding remarks 

In this paper, we have proposed a new type of reduction involving an antisymmetric 
constant matrix; this reduction relates one matrix variable with another in a system of 
coupled matrix PDEs. In the particular case of vector variables, it enables us to obtain the 
integrable vector PDEs having Sp(m) as their symmetry group. Our approach has been 
proven to be applicable to both continuous and discrete systems. For the particularly 
interesting systems such as ( j277j) . ( 12^221) . (PX341) . ( 1377|) . ( 13T201) . and ( 13T251) . the one-soliton 
solutions are derived from the (discrete) linear integral equations of the Gel'fand-Levitan- 
Marchenko type. The solutions are clearly invariant up to a redefinition of the soliton 
parameters under the same symmetry group Sp(m) as that of the original systems. The 
behavior of the solitons reflects an interesting characteristic of these systems; the total 
"particle number" of the system, e.g., ||it|| 2 da; in the continuous case, is, in general, 
not conserved and varies in time. As a result, each soliton exhibits an overall vectorial 
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oscillation, in addition to an internal oscillation among the vector components. A de- 
tailed investigation of the multi-soliton solutions would be an interesting and promising 
path toward the construction of a new class of set-theoretical solutions with symplectic 
invariance to the quantum Yang-Baxter equation (cf. ref. HSJ). 

One natural question arising from the results of subsections 12.11 and 12.21 concerns the 
integrability of a system of the following general form: 



dui d 3 Ui 
dt dx 3 



Yl C i k ( 



duj duk 
-^-Uk ~ uj-^— 
ox ox 



dx 



u k -u 

\ ox 

l<j<k<M 





Qui 
dx 



E fduj du k \ 
jk \dx Uk ~ Uj ~dx) 



Li: 



X<j<k<M 



1,2,..., M, (4.1) 



where a and b are constants. Using the Mathematica package "InvariantsSymmetries.m" [46J, 
we searched for the cases wherein system ( 14. ip with M = 2 can possess higher polynomial 
conservation laws and/or symmetries of the prescribed orders. The result was null in that 
within the limitations of our computer's memory and CPU performance, we could detect 
no integrable case except the already found two cases a = and 6 = 0. Our future task is 
to rigorously prove (or disprove) the nonintegrability of system (14.11) in the case ab ^ 0. 

We would like to explain how the semi-discrete systems of the form d t Un + — 
^n-i) + {nonlinear terms} = presented in section [3] can be related to the third-order 
PDEs of the form dj<Ui + dj^Wj + {nonlinear terms} = in a continuous limit. In fact, the 

asymptotic expansion with respect to the space interval A, d t Un +(u^ +1 — w^-i) — d t u^ + 
2Ad x u^ + ^A 3 d x , u^ + 0(A 5 ), implies that a Galilean plus scaling transformation such as 

d t + 2Ad x =: \A 3 d T , d x =: d x , ~ O(A^), or equivalently, u^(t) ~ A\u, (A(n - 2t), AH/3) 
must be performed. This is a commonly accepted technique (see, e.g., refs. HUH?]), an d 
the solutions of such a semi-discrete system generally have the same structures as those 
of the corresponding continuous system. However, one can also obtain further "natu- 
ral" space discretizations that directly arrive at the third-order PDEs in the continuum 
limit, without resorting to the Galilean transformations. This is achieved by considering 
a proper linear combination of the original semi-discrete system and a higher symmetry 
of it [17JHH]. Note that this "improvement" results only in a minor change in the time 
dependence of certain parameters in the solutions, while a semi-discrete system obtained 
in this manner usually appears rather complicated and less attractive than the original 
one. Therefore, we do not pursue such a direction in this paper. 

We have concentrated on the reductions of the third-order (cu cx k 3 ) flows of the deriva- 
tive NLS (DNLS)-type hierarchies as well as their first negative (u> oc k' 1 ) flows, in both 
the continuous and semi-discrete cases. The feasibility of the reduction is based on the 
fact that the cubic terms in the evolution equation for q and those for r in such sys- 
tems, e.g., (12. 2p or (I2.19p . have opposite signs. This is in contrast with the case of the 
corresponding flows of the matrix NLS hierarchy, e.g., the (non- reduced) matrix com- 
plex mKdV equation (12.101) that permits various reductions, including R = A\Q T A 2 with 
A\ = Ai and A \ = A 2 or Af = —A\ and A\ = —A 2 , but not the reduction R = CQ T or 
R = Q T C with C T = —C. However, it should be noted that the matrix DNLS flows are 
not the only class of systems that permit the reductions of the latter type exploited in 
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this paper. As an illustrative example, let us consider the matrix generalization of the 
Yajima-Oikawa system [49J (cf. refs. 1501455]) . 



iQ t2 + Q xx -PQ = 0, 
iR t2 - R xx + RP = 0, 
iP t2 + 2(QR) X = O, 

and its third-order symmetry. System (14. 2 j) possesses a Lax pair of the form 



(4.2) 







' -h 




" 


Q 


p ' 










u = 


K 





+ 








R 














h 




_ h 




















' 




' 


Q 







" 


Q x 


QR 


V = 


K 2 


h 


+ C 








-R 


+ i 


R 






























Q 






(4.3a) 



(4.3b) 



This implies that the substitution of (14. 3 j) in the zero-curvature condition (12. 4 p yields 
( 14.21) . The matrix Yajima-Oikawa system (14. 2 1) allows the standard reduction R = BQ\ 
= P, = —B, B t = B x = O. The third-order symmetry of (14. 2 p reads as 



Qt 3 + 4Q XXX - 3P X Q - 6PQ X - 6QRQ = O, 



Rt 3 + 



3RP X - 6R X P + 6RQR = O, 

= 0, 



(4.4) 



Pt 3 + P xxx ~ 3(P 2 ), + 6(Q X R - QR X 
and the corresponding Lax pair is given by (I4.3ap and 

Ah 



v = iC 



+ 



o 



O 4Q 4P 
+ C 2 O O AR +iC 
4/x J [Ah O O 

P X + 2QR -4Q XX + 2PQ -P xx + 2P 2 + 2{QR X 
-ARQ -AR XX + 2PP 



2P 4Q X 2P X 

O O AR X 

O O 2P 

— Q X R) 



AR X 
2P 



-AQ 3 



-P x + 2QR 



System (I4.4p allows an extension of the typical reduction considered in this paper, that 
is, R = CQ T and P T = P or R = Q T C and P T C = CP, where C is an antisymmetric 
constant matrix. In particular, the vector reduction in the former case changes the matrix 
system (I4.4p into the following system [5T1I56] : 



du i , d 3 Ui 
— - + 4 

dt dx 3 

dp d 3 p 
dt dx 3 



dp 
ox 



dui 
dx 



= 0, 



1,2 M, 



l<j<k<M 



d 2 u-i 
dx'- 



Un 



d 2 u k 
dx 2 



0. 



It is noted that this system is a modification of the triangular system comprising the KdV 
equation and time part of the associated linear problem due to the addition of the last 
summation term. 
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